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Detailed physical processes of magnetic field generation from density fluctuations in the pre- 
recombination era are studied. Solving Maxwell equations and the generalized Ohm's law, the 
evolutions of the net charge density, the electric current and the electromagnetic field are solved. 
Unlike most of previous works, we treat electrons and photons as separate components under the 
' assumption of tight coupling. We find that generation of the magnetic field due to density fiuctu- 

, ations takes place only from the second order of both perturbation theory and the tight coupling 

approximation. 
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I. INTRODUCTION 



fSl ' Magnetic field generation from density fluctuations in the pre-recombination era (T > 0.3 eV, where T is 

Q ai, " " " " " " 



temperature of the universe) has been investigated intensively by many authors [ll, [2, M 0, S H, S B 111- 
Density fluctuations were generated quantum mechanically during the inflationary epoch, evolved linearly during 
the radiation dominated era, and then acted as the seed for the anisotropy of cosmic microwave background 
1 ■ and large scale structure of the universe. Three important ingredients beside dark matter and neutrinos are 
O ' photons, protons and electrons They basically behaved as a single fluid until the recombination epoch due to 
-j-J , strong coupling through Thomson and Coulomb scatterings. However, because photons preferentially interact 
^ ' with electrons rather than protons, there must be tiny but finite deviation of motion between electrons and 
I— i] protons, that is, the net charge density and the electric current. These are what generate magnetic fields before 
recombination. 

\ This mechanism has attracted considerable attentions because it could give the seed fields for galactic magnetic 
t fields. Galaxies are observationally known to have magnetic fields of order 1 /xG while the origin has been a 
^vq ■ great mystery in modern astrophysics. It is usually considered that if a galaxy has tiny "seed" magnetic fields 
\^ : at its early stage, various hydrodynamical and/or plasma instabilities would amplify the seed fields. This is 
' known as the dynamo mechanism. Accordingly, the problem to find out the origin of galactic magnetic fields 
, reduces to that of seed magnetic fields. Although there have been many mechanisms proposed on the origin 
' of seed fields, the scenario considered here, which is magnetic field generation from density fluctuations, has a 
f — " great advantage compared to other scenarios that it can give a robust evaluation of generated magnetic flelds. 

This is because density fluctuations themselves have already been measured accurately and a theoretical tool 
^ \ to discuss them, cosmological perturbation theory, has been established firmly. Note that second order density 
•rH , perturbations, which are the next order of linear perturbations, are necessary to be considered for generation 
/\ ' of magnetic fields. It is rather lengthy and complicated to solve second order density perturbations and this 
^ , is why reliable quantitative estimation of seed magnetic fields has not been appeared until recent. For more 
- - ' information on seed magnetic fields and other scenarios, see a comprehensive review [lol |. 

So far, the amplitudes of the generated magnetic fields estimated by different authors are roughly consistent on 
the horizon scale at recombination, although they differ by several orders on smaller scales [J S 0, S S 0, @] • 
Moreover, the previous studies have focused exclusively on the magnetic field but not on other electromagnetic 
properties such as the electric field, the net charge density and the electric current. In order to understand the 
physical processes of magnetic field generation, however, it is undoubtedly important to consider all of them 
consistently. For example, one may (wrongly) conclude that magnetic fields cannot be generated because tiny 
deviation of motion between electrons and protons, which induce magnetic fields, could not be maintained due to 
the shorter timescale of Coulomb scattering between protons and electrons than the one of Thomson scattering 
between photons and charged particles. One will find out why this argument is not true later in this paper. 
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The main purpose of this paper is to clarify all physical processes working on generation of magnetic fields 
via density fluctuations and to make a physical interpretation of the results of our previous papers [1, S H, [§1 • 
We will solve Maxwell equations and the generalized Ohm's law to express electromagnetic quantities in terms 
of Thomson scattering term, which is an external force from photons and proportional to velocity difference 
between photons and charged particles. 

As we mentioned earlier, photons and charged particles basically behave as a single fluid such that density 
fluctuations of them evolve together. Because no magnetic field generation takes place in a single fluid limit, 
we need some deviation from that limit. But how large deviation do we need? To answer this question, it is 
helpful to consider the tight coupling approximation This approximation is based on the fact that the 

interaction timescale, Tint, between two fluids is much less than the dynamical timescale, Tdyn- The coupling 
between two fluids is tighter for a smaller interaction timescale and in the limit of Tint/Tdyn 0, we have exact 
tight coupling and no velocity difference between two fluids. Therefore it is useful to expand physical quantities 
with respect to the tight coupling parameter Tint /"Tdyn in order to estimate deviation from a single fluid. 

Let us give a more specific argument. The scattering timescales for Thomson between photons and charged 
particles, and Coulomb scatterings between protons and electrons are, 

rT = ^ -^2x10^ sec ('-±iy\ (1) 



■me . ,„_3 (l + z 



^- = ^-4^i«"'^^^llo^j '''' 

where nUp is the proton mass, ctt is the Thomson cross section, is the photon energy density, rUe is the 
electron mass, e is the electric charge, rie is the electron number density and z is the redshift. Here -q is electric 
resistivity as, 

l + /In A 



where In A is the Coulomb logarithm. 

On the other hand, the dynamical timescale of fluid may be thought as that of acoustic oscillations. For a 
given Fourier mode whose wave number is fc, the timescale of acoustic oscillations can be written as Tac — l/kcg, 
where Cg is the sound velocity. In the radiation dominated epoch, Cg = c/\/3, where c is speed of light which 
we take unity hereafter. Accordingly the timescale Tac can be approximately written as k"^. This simply 
means that a Fourier mode begins to acoustically oscillate roughly when it enters the horizon. This acoustic 
oscillation remains until the fluctuation of the scale damps away due to the finite mean free path of photons. 
This phenomenon, known as Silk damping, occurs at the scale of the diffusion length of photons, kg-^^. Deviation 
from tight coupling also occurs due to the finite mean free path so that magnetic field generation will be effective 
at around the Silk scale, rather than the horizon scale. Thus, the dynamical timescale can be thought to be 
roughly the Silk scale, 



l + z 

f05 



-5/2 



fcsul-Je^-8xfO«sec(^) , (4) 



aTTi 

where Tcos is cosmological time scale, which is the inverse of Hubble parameter, H, 



Tcos = « 5 X fO^ sec ( ^) . (5) 



Then we can have a small parameter for each scattering as, 

-3/2 



Tdyn V / Tdyn V 

where we evaluated the dynamical timescale as kg^l^. The former is the tight coupling parameter for Thomson 
scattering and determines the magnitude of velocity difference between photons and charged particles. On the 
other hand, the latter is for Coulomb scattering and related to the magnitude of the net charge density and the 
electric current. However, as we will see later, deviation of motion between protons and electrons is suppressed 
more than expected from the tight coupling approximation because protons and electrons are coupled not only 
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by Coulomb scattering but also by the electric field. Further, it turns out that deviation of motion between 
protons and electrons is always much smaller than that between photons and charged particles, even though tt 
becomes smaller than tq for T > 10 keV. 

Finally let us give two more important timescales. One is the inverse of plasma frequency, 

and the other is the magnetic diffusion timescale, 

= % -7x10- sec (Iti)"''. (8) 
We see the following hierarchy for various timescales: 

TT, rc, 77, UJp^ < Tdyn < Tcos < Tdiff . (9) 

This hierarchy tells us that we can neglect the diffusion of the magnetic field due to electric resistivity. Also, 
because we focus on the dynamics of the scales much less than the horizon scale, general relativistic effects are 
expected to be unimportant. Therefore Newtonian treatment will be sufficient and we neglect the cosmological 
expansion as well. This treatment will make our analysis rather qualitative but clear, and is appropriate for 
our purpose. 

This paper is organized as follows. In the next section, we derive the generalized Ohm's law and an equation 
for the velocity difference between photons and charged particles. Then in section lllll combining the Ohm's 
law with Maxwell equations, we solve electromagnetic quantities under the assumption of the Thomson term 
regarded as an external force. The evolution of the Thomson term is investigated in section ITVl using the tight 
coupling approximation for Thomson scattering, and combining with the result of section [Till electromagnetic 
quantities are expressed by conventional quantities such as the photon density fluctuation. Finally we give 
discussion and summary in section fVl 



II. EQUATIONS OF MOTION 



We start from Newtonian equations of motion for fluid densities of photons, protons and electrons, neglecting 
the cosmological expansion and the pressure of charged particles. 



4 1 

-p^ [dtv^ + {v^ • V) Ey] = - - [Vp^ + V • (p^n^) 



rUpUp [dtVp + [vp ■ V) Vp] = enp{E + Vp x B) + (7^^^ + - mpUpV^, 
meUe [dtve + {ve ■ V) V,] = -en,{E + VeX B) + Ci'^^ + Ci^'> - me?leV$ 



(10) 

(11) 

(12) 



where V is a derivative with respect to spatial coordinate, Va (a = 7,p, e) are fluid velocities, Il-y is anisotropic 
stress tensor of photons, rip and rif. are proton and electron number densities, respectively, E and B are electric 



and magnetic fields, respectively, and $ is gravitational potential. Here, C^^^'' and are collision terms for 



=;(c) 



Thomson and Coulomb scatterings between i and j particles, respectively 



which are written as 



arfipPj 



{vj - Vp) - -Vp ■ 



de) ^ -ciV = - 



^7e 



£7 



-arriePy 



{V-y 



1 



<5p? = -Cip^ = -e^npUeTjiVp - Ve). 



(13) 

(14) 
(15) 



We will rewrite these equations of motion in terms of center-of-mass and relative quantities of charged particles, 
defined as. 



rib 



Vb 



ip 



■ fine 



npVp + PUeVe 
Up + PUe 



SVr, 



(16) 
(17) 
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and conversely, 



Vb 



rie = Tib 

f3 SUpe 

~ (l+/3)2 Hb 



1 - 
SVv 



' P 



Vb 



1 



(18) 

(19) 



where /? = m^/inrip. In terms of the new variables, the net electric charge density and the electric current can 
be written as. 



p = e{np - Tie) = eSupe, 
j = e{npVp - UeVe) = e 



UbSVpe + SUpeVb 



1-/3 
1 + /3 



6npeSVpe 



rib 



(20) 
(21) 



Then, let us rewrite the equations of motion, Eqs. (fTU|) - (fT^ . in terms of the center-of-mass and relative 
quantities, Eqs. p6|) and (fT7|) . and Sv^b = — Vb- We will keep only linear terms in Supe and Svpe, while we 
will keep all nonlinear terms in Svjb- The neglection of higher order terms in Supe and Svpe will be justified 
later when we solve all equations and find Supe/rib and \6vpe\ are much smaller than |(Su-yb| for temperatures of 
interest (rrig > 0.3 eV). From merieX Eq. pip —mpUpX Eq. p2|l . we obtain an equation for 6vpe, 



e(l + /3) 



dtSVpe + {vb ■ V) SVpe + [SVpe ' V) Vb 

1 + Z?"* CrTP7 



E + Vb X B — I ^n ^Vpe X i? 



1+/? 



eubTi 



(l + /5)2 e 



1 - gTP7 

1 + /3 e 



(JiXyf, - -Vb ■ 



(22) 



which can be regarded as the generalized Ohm's law. Equation of motion for baryon fluid is obtained from UgX 
Eq. (dH) +npX Eq. ^ 



dtVb + {vb ■ V) Vb 



(l + /3)mp 



5vp(, X B + 



1 + P Trip 



(1 + 0^) {Sv^b-]vb-n^] + 



^Svp, 



On the other hand, Eq. (fTO|) can be rewritten as 



dtV.y + (W-y • V) t/-y 



1 Vp^ + V- (p-yH^,) 3 



4 

and from Eq. ([24]) - Eq. we obtain. 



(1 + P^) [ 6v^b - -^Vb ■ 



1-/3%. ■ 

TTp^''^\ 



V$. (23) 



V$, (24) 



dtSv^b + [v-y ■ V) Sv^b + {Sv-yb • V) — {5v^b ■ V) 



1 Vp-y + V ■ (p^n-y) 

4 p^ 
1 + i? crTP7 



1 + /3 TOp 



(i + /?2) Uz/^,--«,.n^ + 



i-p^ 
i + p 



SVr, 



(25) 



where, 



_ 3(mp + me)nb 

^ = 'a 

4P7 



4 X 10' 



l + z 

105 



(26) 



Among these equations, Eqs. (P^ . and ([^5]) can be chosen as independent equations which describe the 
motion of the three fluids. In the context of evolution of cosmic microwave background (CMB) anisotropics, only 
Eqs. and ([^5]) without Svp^ terms have conventionally been considered. This is based on the assumption that 
protons and electrons are tightly coupled through Coulomb interaction. As we will show later, this assumption 
is valid in the sense that the Svpe terms in Eqs. (|24p and are negligible compared to other terms. However, 
it is obviously impossible to argue electromagnetic properties of the early universe in this approach, i.e., the 
complete tight coupling limit or taking only the zeroth order of the tight coupling parameter. 

In the conventional approach, all physical quantities are expanded according to cosmological perturbation 
theory. At the zeroth order, the universe is homogeneous and isotropic so that all the vector quantities vanish. 
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Further, because we are neglecting cosmological expansion densities and resistivity are constant both in time 
and spatial coordinates: = p^^\ Up = Ue = ni, = n'^^ and 77 = 77^"^. At the first order, deviations 
from homogeneity and isotropy are taken into account and typical magnitude of the deviations is about 10~^. 
Density and tensor perturbations have nonzero values which depend on positions as well as time if they were 
once generated during the inflation era. Vector-type perturbations, namely divergenceless vectors such as the 
magnetic field and fluid vorticities, are absent even at the first order since solutions of their perturbations are 
only decaying mode. They can exist only if we consider the second order. Accordingly, we have, 

{t,x)^ + (t , + (t, af) + • • • , (27) 

m{t, x) = nf^ + n^^^ {t, x) + ) {t,x) + -- - , (28) 

T]{t, x) = T]^°^ + 77(1) {t, x) + r;(2) [t,x) + -- - , (29) 

i?(t,f) = i3(2)(t,f) + --- , (30) 

and other quantities start from the first order, although fluid vorticities are absent at the first order, i.e., 
V X v^'' — 0. In this article we will consider up to the second order in cosmological perturbation. We see in 
the above equations of motion that the Lorentz force term, Vh x B, and the Hall term, Svpe x B, are neglected 
because they only appear from the third order. 

III. SOLVING MAXWELL + OHM 

In this section, we solve Maxwell equations and the generalized Ohm's law to obtain electromagnetic quanti- 
ties, p,j, E and B. Maxwell equations and the charge conservation law are written as, 

y -E^ eSupe, (31) 
dtE = V X S — e{n(,Svpe + SupeVb), (32) 
dtB ^ -V X E, (33) 

dtSUpe + V • {nbSVpe + SUpeVb) = 0, (34) 

and the generalized Ohm's law is obtained from Eq. ([22]) as, 

E = 2 a\ [dtSVpe + {Vb ■ V) SVpe + {SVpe ■ V) Vb] + eUblJcffSVpe + C. (35) 



Here ?7eff is an effective electric resistivity which includes contributions from Thomson scattering. 



1 + (TtPj 



1 



(1 + (3y (}tt 



(36) 



Note that the second term in Eq. is dominant for T > 100 cV. On the other hand, C is the Thomson 

scattering term which is regarded as an external force in this section. 



C.'-^^-f^{sv,b~\vb-U,]. (37) 



The purpose of this section is to express p,j,E and B in terms of C up to the second order in cosmological 
perturbation. It is convenient to decompose C into scalar and vector parts, 

C^Cs + Cv, (38) 

where the scalar part Cs can be written by a gradient of a function and Cv is divergenceless, V • Cy = 0. As 
we stated above, Cy is absent at the first order in cosmological perturbation. 
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A. First order 



At the first order, the magnetic field is absent and the equations are, 



at<5nW + nl*'V • = 0, 



■pe I ' "b 
Tor 



pe 



e{l + f3) 

Let us compare the first and second terms in r.h.s. of Eq. (|42p . 



niedtS^fp^ /e{l + P)\ 



krc - 5 X 10-- (^) 
^/3rT - 10-- (^) 



(l + z< 10*^) 
(l + z> 10^) 



(39) 
(40) 
(41) 
(42) 



(43) 



where we evaluated the time derivative by the wavenumber. Thus, we can neglect the time derivative term 
when we consider cosmological scales at temperatures T < m,,'. 



(44) 



Actually, this is the leading-order tight coupling approximation for Coulomb scattering. Taking the divergence 
of this equation, we have. 



'/off ^f'-pe 

We can neglect the first term of l.h.s. because 



\5n^ph 



- 10-" (i^) m) 

-4xl0-i«(^)(i±^) 



7/2 



(l + z< 10**) 
(l + z> 10**) 



(45) 



(46) 



thus we have, 



^ e 

Substituting this into Eq. (|41]) . we can solve for the velocity difference, 



(47) 



(48) 



In general, the rotational part of 6vl]) cannot be determined from Eq. (|^T|) . but it should vanish at the first 
order. Then, we obtain the electric field from the Ohm's law as. 



where we neglected ly^^dtC'^^^ term. Thus, at the first order, we have. 



pe 



(49) 

(50) 
(51) 
(52) 



and we see they also satisfy Eq. Correspondingly, the net charge density and the electric current can be 

expressed as, 



(53) 
(54) 
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Before we proceed to the second order, let us explain the meaning of the approximations in Eqs. (j43p and 
If we take the divergence of the Ohm's law ((1^ . without neglecting the time derivative term, we have, 

(55) 



(l+/3)^2 



Eq. (j55p describes the dynamics of charge separation, 6n^J , and can be seen as an equation of damped oscillation 
with an external force. There are two key timescales. One is the timescale of oscillation, lu~^, and another is the 
damping timescale, l/{LUpr]cs) tq. Even though the charge separation and its time derivative are zero initially 
(outside the horizon), the source term induces the oscillation whose timescale is cOp^. Then the oscillation is 
damped within the timescale rc , and charge separation relaxes into the equilibrium value which is nonzero due 
to the presence of the source term, eSn^pJ — V • C^^^. Eq. ((5^ tells us that the force from photons balances 
with electric field in this state. Because we are focusing on the dynamics of cosmological timescale, it is enough 
to consider the equilibrium state. The first approximation, Eq. (f43|) . corresponds to the neglection of the first 
term in l.h.s. of Eq. (|55p . that is, the neglection of plasma oscillation. Likewise, the second approximation, Eq. 
(Uni, corresponds to the neglection of the second term in l.h.s. of Eq. Note that this neglection leads to 

the absence of magnetic diffusion, that is, magnetic diffusion is important only for very small scales. Thus, both 
approximations are valid in our context. Similar approximations will also be used at the second order below. 



Second order 



At the second order, the equations are. 



9*^(2) = -Vx ^(2)^ 

m = 



e(l + (3) 
io) , 

b '/eft 



dtSi'J + {i'^ ■ v) s4l^ + (5.-<i) 



"6 '7cff / 



As in the case of the first order, the terms in the bracket [• • • ] of Eq. (|60p are suppressed by the factor 
and can be neglected. Further, the last two terms in the parenthesis (• • • ) of Eq. (jHO)) can be ne; 
compared to C'^' because. 



(1) 

(0) (0)^c-<l) 
'/off (0)">e 



„(1) 

_(0) (0)^.-^1) 
^eff 



^(2 



thus we have simplified second-order Ohm's law as, 

As we did at the first order, we take the divergence of this equation, 

and neglecting the first term, we obtain, 



Then we can solve for SvpJ from Eq. ([59]) as 



1 



.(1) 



(56) 
(57) 
(58) 
(59) 



(60) 

in Eq. 
elected 



(61) 

(62) 
(63) 

(64) 
(65) 
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where D^"^"^ is an undetermined vector. From Eq. (j60p . we have, 

^(2) =(7(2)+en["'77(0)Vxi3(2). 

Here we neglected terms of order ki]^ . The magnetic field can be obtained from Eq. ^ 

= - fdtVx 



dt 



(66) 

(67) 
(68) 



and the vector Z)^^) determined by Eq. ([S7|) : 



,(0) 



V X z3(2) 



dt en[°^?7^°V x V x V x D^"^^ 



= - / V X V X „ 9^(7(2)^ 



(69) 



Evaluating time and spatial derivatives as wavenumber k and time integration as 1/fc, we pick up only dominant 
terms: 



V X z3(2) = - 



1 



,(0) 



V X V X C(2) 



1 n^(2) 



Thus we obtained all the second-order quantities in terms of C*^^)^ 

1. 



5n(2) 



5w(2) 



1 



en 



(0) 



llc(i)') + (V • + / dt V X V X 



i?(2) = 



i3(2) VxC(2), 
and correspondingly, 



■(2) = -ac(2) - i dtv xv x& 



(70) 

(71) 
(72) 

(73) 
(74) 

(75) 
(76) 

(77) 
(78) 

(79) 
(80) 

up to the second order. This is one of our main results. Here it should be noted that the vector-type perturbation 
is absent at the first order so that V x C*^^^ =0. 

Now we expressed the electromagnetic quantities in terms of the Thomson term C. We see that the electric 
current has two contributions in Eq. (|78p and the first and second terms balance with the displacement current 
and rotation of the magnetic field in (|57p . respectively. Since the Thomson term C is the velocity difference 
between photon and baryon fluctuations, it suffers from Silk damping when k'^^^ exceeds that is, C — >■ 
for < fc^iJi^. The above solution tells us that the electric current and the magnetic field keep non zero values 



Combining the first and second order results, we have, 

j^^dtC ~ I dt\/ xS/ X C, 
B = - I dtW xC, 
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even for C — while the net charge density and the electric field vanish. This is because the diffusion timescale 
of the magnetic field is so much larger than the dynamical timescale that the magnetic field, once generated, 
would not damp away and be maintained by the residual electric current (the second term in Eq. ([75])) even 
after the source term would disappear. 

Before concluding this section let us comment on the earlier work of 7']. The difference in the treatment 
between theirs and ours is that we have included the displacement current in Eq. (|57p while they have omitted 
it. We found above that the displacement current leads to the second term in ()78p and it is what maintains 
the magnetic fields after their generation. We attribute the difference in the resultant magnetic fields at small 
scales between [7j and ours to this contribution from the displacement current. 



IV. EVOLUTION OF THOMSON TERM 



In this section, we follow the evolution of the Thomson term to discuss magnetic field generation in the 
language of the tight coupling approximation. This approach was first studied in Q and it was shown that 
magnetic field generation is absent at the first order in the tight coupling approximation and starts from 
the second order. We will confirm this fact and, furthermore, solve the equations explicitly to express the 
electromagnetic quantities in terms of the conventional quantities such as the density fluctuations of photons. 
For simplicity, we will ignore the anisotropic stress of photons hereafter. This treatment makes our analysis 
rather qualitative but it would be helpful in order to see the dynamics of the tight coupling approximation 
clearly and interpret the numerical analysis of our previous papers @, |8| . 

When we consider magnetic field generation, a particularly important quantity is the rotation of the Thomson 
term, 



l + (3 



L Pi 



(81) 



Here we ignored the anisotropic stress of photons. Thus, we need to solve for Sv-yb up to the second order 
in cosmological perturbation. In Eq. ([81]) the first term of the r.h.s. is the vector product of the gradient 
of the photon energy density and the velocity difference between photons and baryons. This was called the 
"slip term" in IJ and evaluated numerically. On the other hand, the second term is the vorticity difference 
between photons and baryons and has not been evaluated so far. Although it is hard to compute this term in 
the framework of general relativistic perturbation theory, we can evaluate it in our simplified formalism. 
The evolution equation of 6vjb, Eq- (ESI), without the anisotropic stress of photons is. 



dtSvjb + {v-y ■ V) Svyb + {Sv-yb • V) tty — {Svyb ■ V) Sv^b = — 7 



1 Vp^ I + R arp-y 



4 1 + P nip 

We can compare the relative importance of the collision terms from the result of the previous section. 



{l+P^)dv^b + \ , Svp 

(82) 



1 + 13 



\6vp^\ _^ aTP,k _ k ^^^^^^-.s(Jl_)(1+^V^\ (83) 



\Svyb\ e^Ub Lu'^TT \ksiikj V 

Therefore, we can neglect the Svpe term and we have, 

dtSv^b + iv-V) Sv^b + {Sv^b ■y)v- {Sv^b ■ V) 6vyb = -7— - 1^(1 + R)^^Sv^b- (84) 

We have also confirmed the conventional assumption in the context of evolution of CMB anisotropics that 
protons and electrons are so tightly coupled that we can treat them as a single fluid as long as we consider the 
dynamics of photons and baryons. 

We will solve Eq. ([84|) by employing the tight coupling approximation [ll|. As we mentioned in the 
introduction, this approximation makes use of the fact that the scattering timescale, tt, is much shorter than 
the dynamical timescale, k~^, so that the deviation of the motion of the two fluids is very small. The expansion 
parameter is the ratio of two timescales. 
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Here we define the deviation of photon and baryon distributions from adiabatic distribution by, 

Pj=Pj{l + A^), nfc = nb(l + Afc), (86) 

respectively. These two quantities, and Af,, are assumed to be small and we expand them by the tight 
coupling parameter, Eq. (|85p . as, 

A^ = A(^) + A(^^) + • • • , A, = a[') + a['') + • • • . (87) 
Fluid velocities are also expanded as, 

^^ = ^+^W+viii) + ... , + +4"^ + ■■■ , Sv^t = Sv^^i^ + Sv^^i'^ + ■ ■ ■ , (88) 

where v is the common velocity of photons and baryons at the zeroth order and Sv^^^ — v^^"^ ~''^\ (* = 1 1 Hi ' ' ' )■ 
It should be noted that this expansion is independent of cosmological perturbation and each term in Eqs. ([57)) 
and (|88p can be expanded with respect to cosmological perturbation, e.g. a(^) = a(^'^) + a(^'^) + .... 

Our purpose here is to express dv-yb and then the electromagnetic quantities explicitly in terms of the con- 
ventional quantities such as ^7^^. 



A. Continuity Equations 

To solve the equation of motion ((M)) . we need to consider the continuity equations of photons and baryons, 

dtp^ + {v^ ■\7)pj + -v^^O, (89) 

dtUb + {vb • V) nf, + nfeV -Vb^O. (90) 
At the zeroth order of the tight coupling approximation, these equations reduce to, 

dtp^ + {v-V)p^ + ^^\/ ■v = 0, (91) 

dtfib + {v-V)nb + UbV • ?7 = 0, (92) 

which are combined to obtain, 

{dt+v-^)l^^)^^Q. (93) 

This shows that fluctuations of photons and baryons behave adiabatically at this order, as we expected. 

As we will see later, we need the relation between A-y and Af, at the first order both in the tight coupling 
approximation and cosmological perturbation to discuss magnetic field generation. At this order, the continuity 
equations for photons and baryons are. 



a,A(^-i) + |V • vil^^^ = 0, dt^i'''^ + V • 4^'^) = 0, (94) 
respectively, and then wc have. 



Ar = + / V . (95) 



B. Equation of Motion 

Substituting the expansion Eq. (j87p into Eq. (|84p . we have, 

dt6v^b + {V'j ■ V) 6v^b + {Sv'jb • V) zXy — {Sv^b ■ V) Sv^b — (^—r-^ + VA^ j — i^Sv^b, (96) 
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where v is the collision frequency, 

V = 

1 + /3 



(l + i?)^^^(l + A^) 



1 



R 1 + Ab 



1 + R 1 + i? 1 + A, 



(97) 



and barred quantities are the zeroth order of the tight coupling approximation: 



114(1+ i?)^^ 



R: 



3(TOp + me)n}, 



At the zeroth order in tight coupling approximation, the equation of motion reduces to 



(98) 



= 



IVp^, _ (/) 

vbv 



4 P. 

This equation can be solved up to the second order in cosmological perturbation: 



(99) 



6v, 



{1,1) 

76 



1 vp-w 



(100) 
(101) 



We can show that the rotation of the Thomson term ([5T|l . which is the source term of the magnetic field, 
vanishes at this order, denoting the following relation obtained from the adiabaticity condition, Eq. (|93p . 



m ~ 4(l + i?(o)) pW' 



(102) 



Let us now consider the next order of the tight coupling approximation in order to argue magnetic field 
generation. The first order equation is. 



(103) 



and this can be solved as, 



1 



1 dt^p' 



-rb 4P(0) ' 4(P(0))2 pio) 



(1) 

7 



(104) 



6v 



{11,2) 

76 



1 



4P(0) 
1 



4(j>(0))2 



77(1) 

(1) 



d^P_ 



P7 



P7 



. Pf ^ 



Pt^ 



• V U 



(1) 



dtp\ 
. pf 



Pt^ 



(105) 



where we substituted Eqs. (|100p and poip . We see that we have non-zero contribution for the slip term and 
vorticity difference at this order as 



P7 



(1) 



Vp. 



Pi 



4(i?(0))2 ^(0) 



X 



(1) 

7 



(106) 
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4P(0) 



X a(Ai) 



1 



4(j7(0))2 
1 




4 + 3i?(o) Vp^'^ 



16(P(0))2 l + i?(o) p(0) -(0) ' l6(p(o))2 i + i?(o) -W 

where we used Eq. p02p and the following relation derived from Eq. ([M]) . 



4(l + i?(0)) ^ 



dt 



-(0) 



z^(o) 4(l + i?(0)) ^ 4(l + i?("))i>(0) 
Adding these two contributions, the rotation of the Thomson term can be written as, 



dt 



(107) 



(108) 



V X = 



1 l-/53^Tp(,0) ^(0) 



16(p(o))2 1 + ^ e l + i?(o) pW 



(1) 

■7 



dt 



V(V^P^^^) 



(109) 



Below we show the leading order quantity for electromagnetic quantities. All quantities except the magnetic 
field need just On the other hand, as we saw above, the magnetic field vanishes at the first order in tight 

coupling approximation and the second order terms, i.e., Sif^^'^^ and d'tf^^''^^ are necessary to have a nonzero 
slip term and vorticity difference in Eq. (|8ip . respectively. 



11- (3^ 



1 



nip V^p^"^^ 



4 1 + /32 1 + i?(0) e 
II- 13^ 1 mp dtVp. 



-(0) 

Pi 



(1) 

■7 



4 1 + /32 1 + i?(0) e 
1 1 - 1 



-(0) 

P7 



5(2) = 



4 1 + /32 1 + i?(o) e 
1 (l+/3)(l-/33) 



16 (l + /32)2 



(l + i?(o))3ea^p(«) 
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1 



4 1 + /32 1 + i?(o) e2 
1 1 - ^3 2 nip dtVp. 



-(0) 
P7 



(1) 

7 



4 1 + /32 1 + i?(o) e2„^ 



P7 



-(0) 

L pV 



-(0) 

P7 



2^(1)^ 



dt 



v(v2py 

-(0) 

P7 



(110) 
(111) 
(112) 
(113) 



We see that all quantities are expressed by some background quantities and a single perturbed quantity, p^\ 
In Eq. (|113p . the first term in the bracket is contributed from both the slip term and vorticity difference, while 
the second term is contributed only from the vorticity difference. 



V. DISCUSSION AND SUMMARY 



In this paper, we made a physical interpretation of the results of previous studies on magnetic field generation 
from density fluctuations in the pre-recombination era. This was done by solving Maxwell equations, Ohm's 
law and an equation for velocity difference between photons and baryons. First we expressed electromagnetic 
quantities in terms of the Thomson term and studied their behavior. We saw that timescales for Coulomb 
scattering are so short that charge distribution quickly relaxes into its equilibrium state, which is not charge 
neutrality but charge separation which balances with the external force from photons. It was also shown that 
magnetic field and electric current do not vanish even after the source term disappear in contrast to electric field 
and charge density. Then the Thomson term was obtained by the tight coupling approximation up to second 
order and electromagnetic quantities were expressed by conventional quantities such as the density fluctuations 
of photons. We found that the second order terms in the tight coupling play essential roles for generation of 
the magnetic field. 

Let us give some order-of-magnitude estimation of various quantities obtained in this paper. The deviation 
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of motion between photons and baryons is evaluated as, 



4 ^ 



dv 



3 X 10"^° 



k 



(0) '^^T 



Ik 



1 + z 

105 



-3/2 



10- 



(114) 



This is exactly what is expected from the tight coupling approximation. On the other hand, the magnitudes of 
electromagnetic quantities are, 



en 



(0) 



3 X 10 



4e2nl°) 



'34 



■Silk 



1 fc^ 

^ 1 + .^^ 
105 



10- 



(115) 



This is much smaller than expected from the tight coupling approximation for Coulomb scattering. The reason 
for this can be seen in section lTlI A) where we used two types of approximation. The first is the conventional tight 
coupling approximation for Coulomb scattering, Eq. (|43|) . and the second is Eq. (|46|) . Thus electromagnetic 
quantities are suppressed by two factors. 



X fcryoff 



fc2 



(116) 



which is the suppression factor seen in Eq. (|115p . The second suppression factor can be attributed to the 
existence of coupling due to the electric field. Equations pi4p and (IllSp justify our strategy on the tight 
coupling approximation that we solved up to second and first order for Thomson and Coulomb scattering, 
respectively, dropping the nonlinear terms in (Jn^e and (Supe and keeping the nonlinear terms in 5v^h, and 
Db up to second order. 

It is instructive to know how equations of motion for photons (fT0|) . protons (fTTjl and electrons (fT2|) are 
balanced, each term of the equations at the leading order, 
photons: 



(acceleration) : (pressure) : (Thomson(with protons)) : (Thomson(with electrons)) ^1:1: 0^R : R, (117) 



protons: 

(acceleration) : (electric field) : (Thomson) : (Coulomb) ^ 1 : 1 : : krj, (118) 

electrons: 

(acceleration) : (electric field) : (Thomson) : (Coulomb) ^ (3 : 1 : 1 : krj. (119) 

It can be seen that Coulomb scattering is not important for cosmological scales with hi] <ti 1, which corresponds 
to the fact that magnetic diffusion is absent at these scales. 

In this paper, while we ignored the cosmological expansion and the anisotropic stress of photons, which makes 
this treatment rather not quantitative but qualitative, we could develop a comprehensive treatment of photon, 
electron and proton fluids, and successfully understand how the deviation between these fluids are determined. 
We found out deviation between photons and charged particles is much larger than that of protons and electrons. 
Our analytic treatment is particularly important on scales much smaller than the cosmological horizon where 
we can safely ignore the cosmic expansion and the numerical analysis is rather difficult to carry out. To make 
a precise prediction of the spectrum of the magnetic field, however, we need to evaluate the contribution from 
vorticity difference together with the anisotropic stress of photons, which we will present in future. 
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